THE ALIQUOT CONSTANT 

WIEB BOSMA AND BEN KANE 

Abstract. The average value of log s(n) / n taken over the first N 
even integers is shown to converge to a constant A when N tends to 
infinity; moreover, the value of this constant is approximated and 
proven to be less than 0. Here s(n) sums the divisors of n less than 
n. Thus the geometric mean of s(n)/n, the growth factor of the 
function s, in the long run tends to be less than 1. This could be 
interpreted as probabilistic evidence that aliquot sequences tend 
to remain bounded. 



1. Introduction 

This paper is concerned with the average growth of ahquot sequences. 
An ahquot sequence is a sequence Oq, cti, a2, • • • of positive integers ob- 
tained by iteration of the sum- of- aliquot- divisors function s, which is 
defined for n > 1 by 

s(n) = Y.d- 

d\7L 

d<in 

The ahquot sequence with starting value Cq is then equal to 
ao, ai = s(ao), 02 = s(ai) = s^(ao), 

we will say that the sequence terminates (at 1) if = s^(ao) = 1 
for some k > 0. The sequence cycles (or is said to end in a cycle) 
if s^{ao) = s'(ao) for some k,l with < I < k, where s^{n) = n hj 
definition. 

Note that s is related to the ordinary sum- of- divisors function a, 
with a{n) = J2d\n d, by s(n) = a{n) — n for integers n > 1. 

The main open question in this area can be phrased as: does every 
aliquot cycle remain bounded? That is, does every aliquot sequence 
terminate (at 1) or cycle, or do sequences exist that grow unbounded? 
The conjecture that all sequences remain bounded is often referred to 
as the Catalan-Dickson conjecture. 

The origin of this paper lies in computational work done to test inte- 
ger factorization routines for the computer algebra system Magma |T]. 
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The most efficient known method to compute s{n) uses the multiphca- 
tivity of a and requires the factorization of n. Iterating s provides 
long sequences of more or less random numbers of similar size, and 
this property is useful in testing factorization methods. It was noticed 
that for even starting values the aliquot sequences tend to increase or 
decrease in size fairly slowly, by an amount that seemed constant over 
different starting values, whereas sequences with odd starting values 
usually terminate quickly. 

Around 1996 Andrew Granville [0] furnished a proof for this phenom- 
enon; see Theorem fl3.3p and Theorem fl3.7p below. Further computa- 
tion seemed to suggest that the constant A involved would be smaller 
(but only just!) than 0, but no proof of this was obtained. Recently, 
we were able to obtain estimates that are good enough to prove this 
property. 

Altogether this led to the main result of this paper. 

Theorem 1.1. The geometric mean /i of ^^^^ over all positive integers 
n exists, and equals 

In particular, the aliquot growth factor fi = 0.969 ■ ■ ■ < 1. 

Roughly summarizing, this means: on average, even aliquot sequences 
tend to decrease in size! In some sense this may be taken as probabilis- 
tic evidence in favour of the Catalan-Dickson conjecture. 

This paper is built up as follows. After some preliminaries, we state 
and prove the convergence of the geometric mean for even values; also, 
an expression for the resulting constant A = log /i as a difference of a 
(closely related to the growth of a{n)/n) and /3 is derived. In the next 
section an easy upper bound for a (which involves a sum over all prime 
numbers) is given. The final section is devoted to a lower bound for /3; 
this is trickier, as it involves an infinite sum of terms that themselves 
are infinite products over all primes. 

2. Elementary observations 

Although (J^n) > n for n > 1, all three possibilities s{n) < n, s{n) = n 
and s{n) > n for s do occur: s{p) = 1 for prime numbers, and in 
general s{p^) = 1 + p + ■ ■ ■ + p^~^ < p^ for prime powers; s{P) = P for 
perfect numbers P = 2^~^(2^ — 1) (with 2^ — 1 prime), and s{n) > n 
for n = P -q, where P is perfect and q any odd prime other than 2^ — 1, 
since a{Pq) = 2P(g + 1). 

Besides terminating at 1 after hitting a prime, or ending in a perfect 
number, it is also possible that an aliquot sequence ends in a cycle of 
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length 2 or more: amicable numbers are pairs m, n for which (T(m) = 
m + n = (T{n), hence s{m) = n and s{n) = m, and a 2-cycle is formed. 
Sociable numbers form cycles of larger length (and are known only for 
length 4, 5, 6, 8, 9 and 28; see [10]). For more on these cycles, including 
historical remarks, see also [8]. 
Note that 

= n ■■■+/) 

p^\\n 

is multiplicative, while s{n) = a{n) — n is not. A useful observation is 
that 

(2.1) !^ > 

by multiplicativity of a and since ^ > I for A; > 1. 

Also note that 1 + p + ■ ■ ■ + p'^ (for prime p) is only odd when p is 
odd and k is even. Hence cr{n) — n for odd n will only be even if n is 
a square, and for even n it will only be odd if n is a square or twice a 
square. Hence: unless an accidental square (or twice a square) occurs, 
parity is preserved in aliquot sequences! In this sense s does not behave 
randomly at all. 

In fact, divisibility by (even) perfect numbers also tends to persist, 
and Guy and Selfridge [7] studied other types of persistence as well, 
but we will ignore all but the parity aspect and only consider even and 
odd aliquot sequences separately. 

We are interested in the growth of the sequence n, s{n), s'^{n), . . ., in 
other words, in the question of whether s{n)/n tends to be smaller or 
greater than 1. 

Apparently, Wunderlich (in [15]) was the first to state the following 
result explicitly; the first statement (formulated for a rather than s) 
appears already in [5]. 

Theorem 2.2. Theorem 

hm if^M = !^-l = 0.6449... 

N-*oo N n 6 

1 j:^, s{2n) 57r2 

..mlE^e^ = 1. 0.2337... 

N^oo N 2n-l 24 
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As 

f s{2n - 1) _ ^ 

71 odd 

we find by the same argument usually given for the computation of 
C(2), that 

JV-oo^^iV 2n-l frt x2 4^^ ^ 

a; odd 

from which the whole theorem follows. 

Based on the second statement in Theorem (12.21) . Guy and Selfridge 
seem to have drawn the conclusion that even aliquot sequences will tend 
to grow unbounded (see [7j page 103). Just like a sequence in which 
the terms are alternately multiplied by 5 and by | will remain bounded 
although the average growth factor tends to 2.6, we cannot draw the 
conclusion that even aliquot sequences tend to grow unbounded from 
the fact that the average of s(2n)/2n exceeds 1. What really matters 
is not the arithmetic mean, but rather the geometric mean: 



n = exp f log(s(n)/n) 



\n=l ^ 



n=l 



Remark 2.3. To draw conclusions about the Catalan-Dickson con- 
jecture, one needs more than just the arithmetic or geometric mean 
of a{n)/n. Davenport [3] showed that there exists a continuous func- 
tion of t giving the natural density of t-abundant numbers satisfying 
a{n)/n > t. See [13], and [11], for recent progress on this function. 
This needs to be combined with the persistance of drivers as in [7]. 

A good approximation is known for the value of the function at t = 2, 
implying that s{n)/n exceeds 1 for a little less than a quarter of all n, 
see [1]. 

Finally, note that these questions relate to deep problems, as it is 
known [12] that the Riemann hypothesis is equivalent to the statement 
that a{n)/n is bounded by e''' log log n for all n > 5041; see also [2], [9], 
llU on this connection. 



3. The aliquot constant 



We first show the following result on the geometric mean for the ordi- 
nary sum of divisors function. 
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Proposition 3.1. 



1 y log ^ = A + ol -—] , with A= y a(p) ^ 0.4457, 



where 



(3.2) a{p)= (l--) V — log (1 + - + ■■■ + — ) . 

Proof. Taking the product over all powers p"^ dividing n, with p prime 
and m > 1, we have 



and hence 

n<x n<x p^\n y p'"<x' r\f j 

= a; > — log + O . 

pf^.P'" pijp^^-l) \\ogxl 

For a fixed prime p we have 



, 1 + r = 1 + , 

p[pm _ \^ p[p"^ — 1) 

so 

^ 1 , - 1 V- 1 , 1-1 /P"""^^ 

> — log = > — log 

p^.p"^ p{p"'-l) ptt.p"^ 

P ^ ( p) ~^ ^P P^) ^ ( P^) ~'~ ~'~ ^2;^ 

( J*) (p ^ ( ~'~ p) ~'~ ^ ( P P^) ~'~ ) ~'~ 
and the result follows. □ 

Theorems 13.31 and 13.71 are our main asymptotic results on the growth 
of aliquot sequences. Roughly speaking, they state that the growth 
factor diverges to when considered over all starting values, whereas 
confined to even values it converges, to A. 

Theorem 3.3. 

1 ^ s(n) 

— 5] log ^ = -e~^loglogA^ + 0(logloglogiV). 

^ n=l ^ 
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Proof. As s{n) > 1 for n > 1 we have 

n ~ pi 

for the smallest prime factor pi of n. Thus 

J2 log^^>- logpi(n) > -^logp 1 

l<n<x ^ l<n<x P<x ri<x 

> — ^ logp ■ #{m < — : m has its prime factors > p}. 

p<x P 



By the small sieve: if p = x « then 

[l--|-[l + ^(e— ) + -^ 



#{m<-:g|m^g>p} = FT ( 1 - 1^ - f 1 + 0(e-2«) 
P q<p \ q) P\ 



1 + 



logx^ 

2n \ 1 



plogp y \logp / ) 

Hence 

(3.4) > — e~'''xloglogx + 0(x). 

On the other hand, let M be the set of integers of the form mp < x 
where all prime factors of m are > plogx. For such integers we have 
ui{m) < log(x)/log(plogx), 

a{mp)/mp < (i _ 1 V < (i _ ^ ) i +o ^ ' ^ 



(y{p)IP tz\ qj \ p\ogx) \p\og\ogx^ 



q\m 

SO 



a(mp) 1 / ^ / 1 

1 < - 1 + O ' 



mp p \ Y log log a; 

Then 

log = 2^ log I 

l<n<x 



E log^ = E logf 

<n<x ^ l<n<x \ 
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Let M' be the set of integers mp E M with p < where U is a large 
fixed number. With y = 6^'°^^°^^^ we find: 



logp > ^ logp > ^ logp ■ # < m < — : q\m =^ q > plogx > 

n€M neM' ^ ^ -k [ P J 

y<p<x u 

y<p<xTr 9<Pl°g^- ^ '^/^V \ 6// 

- P ^logx)) 

y<p<x u 

(3.6) > e~'*^x (loglogx + O(logloglogx)) . 

Combining (13.61) and (13. 5p with (13.41) and Proposition (13.11) we obtain 
^ g (fi) 

log = —e~"'x log log X + 0{x log log log x) . 

l<n<x ^ 



□ 



Theorem 3.7. 



^|:iog^ = A+o(i/iogiv) 



where 

( \ 

(2/3,(2) - 1) n i^Ap) 



(3.8) A = a(2)+ «(P)-E- 

p prime i>l 



\ prime / 



with 



(3-9) «(!>) = f 1 - -) E 4 log fi + - + ■ ■ ■ + 4rl 

V pj^iP \ p p J 



as before, and 



l\ ^ 1 A 1 1 ^ ^ 



(3.10) (3jip)= 1--^1 + - + --- + — 

\ Pj ^oP \ P P 



Proof. Suppose that J is sufficiently large. Then 

2n 



fa(2n) \ a(2n) / 



(2n), 



2" fclJ,feV'^(2")/ V(3/2)- 
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using (12. ip . Proceeding as before we get 

(3.11) Elog^-A*x + of^), 

where 

= E + E ^ (log2 + log (i - ^)) • 
We now use the following result. 

Lemma 3.12. Let f be a multiplicative function with < f{p"'') < 1 
with 1 — f{p"^) <^ - for every prime p. Then 



(3.13) 5:/(n)=xn 



n<x 



+ 0((logx 



Proof Let (^(p-) = /(p-) - K i, and ^?(r2) = Dp-^Hn ^^(p'")- Then 

E/w = E n(i+^?(p'")) = 

n<x n<x p"^\\n 

= E E ^?(rf) = E^?(^) E 1 = 

n<X d\n d<X m<x/d 

gcd{d,n/d) — 1 gcd(m,d) — 1 

d<x d<x 

p \ m>i y ^ 



But 



<r.m, \ rri I ^— ^ ^ITl 

m>i y y \ y/ rn>i y 

V Pj^o P"^ ' 

and the Lemma follows. □ 



We apply this to fj{n) = (^)^ then 



fAP" 



P 



m 



3 



a{p"^) 
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for p an odd prime, and 



3-2 



m \ 1 



a{2 



im+1 



But then by the Lemma 



2n 



J my t 



n 



X 



n 



J (3/2)^- V 




+ 0((logx)^) 



X 



-(2/3,(2) 



n /9,(p)| +0((logx) 

\ prime 



Collecting the information we get 

a(2n) ; 



n<x 

and the result follows. 



2n 



log a; 



□ 



Example 3.14. Although most of the rest of this paper is devoted to 
a numerical estimate for A, necessary because of the behaviour of /9, it 
is easy to see the convergence of A. If we sum s{n)/n for the first even 
values of n, up to N , we get the following. 





T.2n<N s{2n) / {2n) 


10^ 


-0.0567457527... 


10^ 


-0.0356519058... 


10^ 


-0.0335201796... 


10^ 


-0.0332873082... 


10^ 


-0.0332626444... 


10^ 


-0.0332598642... 


10^ 


-0.0332595156... 


3.953 ■ 10^ 


-0.0332597045... 



By definition, 

A = 2a(2) + 



p>3 
prime 



4. Computing a 
/ 

( 

V 



E "(p)-E7 



(2/3(2) 



n 

p>3 

prime / 
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In this section we will compute a good approximation and upper bound 
for: 

a = 2a(2)+ ^ a{p). 



p>3 
p prime 



First note that for any prime p 



1-- V — log l + - + ] = y] — log' 



P) n^^lP"^ \ P P"' ) m^lP" \P+---+P 

Since for all m > 1 

p + ■ ■ ■ + p"^ p + ■ ■ ■ + p^ ~ 

and log(l + x) < X for all x > 0, we find that the tail 

is bounded by 

^ '^l + - + 4 + ---VloKfl + ^^^ < p ( ^ ^ 



pAI+l \ p pi I \ I — p _ I \pM+l 

We will denote this bound by 



A{p,M) 



/I \ 2 
p / 1 \ 



1 \pM+^ J 

It also follows that, for any > 1 

E «(p)< E E4iogfi + 4rl<E-iog(i + -)- 

P>iV p>JV rn>lP \ ^' / n>N ^ ^ ^ 

p prime p prime 

But 

V - log f 1 + -) < - log (l + -) dx= r- log (1 + z) dz 

n>N'^ \ n/ Jn X V X/ Jo z 

is clearly bounded by ^. 

Theorem 4.1. For any N > 2 and L, M > 1: 

« < E^log(l + ^ + --- + ^)+2A(2,L) + 

+ E E4i"sf ^!"""!''I l+ E + 

p oddprimc p prime 
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As the sums in this theorem are all finite this gives us an effective 
way to compute an upper bound on a. 

Example 4.2. In the table below we have listed the outcome of some 
computations for a with L = M = 15. These computations were done 
within half an hour (including the primality tests for all primes up to 
10^) on an ordinary PC, using Magma. 





sums 


error bound 


10^ 


0.69^ 


^3072233... 


1.0000093132... ■ 


10' 


-4 


10^ 


0.69^ 


^3162365... 


1.0000931323... ■ 


10' 


-5 


10^ 


0.69^ 


^3169710... 


1.0009313233... • 


10- 


-6 


10^ 


0.69^ 


^3170329... 


1.0093132338... ■ 


10' 


-7 


108 


0.69^ 


^3170383... 


1.0931323384... ■ 


10' 


-8 



Corollary 4.3. a < 0.69831705. 



5. Computing (3 
Our next goal is to compute a lower bound for 



/ 



j>i J 



(5.1) 

where for every prime p 



(2/3,(2) -1) n 



p>2 
prime 



1 I i I J_ J L J_ 



Observe that for j > 1 



P 



m \ J 



1 + E 



1 {p" 



a{p"^) J 



p 



2'm-l\j 



m>l P' 

1-E 



^ (t(p"^)J(t(p'""1)J 
1 Ei=i ii){p'^-^)\p^a{p^~^)y~^' 



m>l 



p" 



a{p"'ya{p'^-^y 
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m>0 

where the multiphcative function z/(n) denotes the number of different 
prime divisors of n, and we set for non-negative m 

( m ^_J_ ( V 

(in particular gj{V) = 1), and for positive m 

j 



while = 1 by definition. 
If we extend our definition to 



I (pcr(p'"-i))'= 



n 



for any positive integer n, we find 

9^ in) = n 9Ap" 



p'" ||n 



where n = Op^HnP™' is the factorization of n. We also define hj for all 
positive integers n by multiplicativity 

h,{n) = n h.ipn. 

||n 

If we define, for composite n, 

(5.2) P,{n) = {-ir^-^g,{n)h,in); 

then for every positive integer n, we find 

n /?.(p) = E/5.H- 

p prime 

This way, the infinite product in the definition of f3, (15.11) . is replaced 
by an infinite sum; note, however, the slight complication caused by 
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the factor 2/3^(2) — 1 rather than f3j{2) in this definition. Since 



(5.3) 2P{2)-1=J:^( \ 



E 9 An 



m>l 



we obtain 
(5.4) 

where 



E 



\ n even 



(5.5) (3*{n) = gj{2^)(3j{no) for n = 2^no, with rio odd. 

To get a fower bound for (3 we will replace these sums by finite 
summations for bounded j and n, and bound the remaining terms. As 
clearly [3j{j>) > for odd p, and so is 2[3j{2) — 1 by (15. 3p . we see from 
the definition that for any J > 1 



/?>E-| E 

j=i J 



nGN 
\ n even 



Define for real e > 



nGN: hAn) > — } . 



Since 



E 

n — l^nQ even 
no>N, no^Sj g 



< 



E 



n 



71— 2^ no even 
no>N, no^Sj , 



n \a[n) 



hjirio) < 



2^nn 



OO 



hj{no) < 



no>N 
odd 



3 J 2Jx=N x^+^ 2eN^ V3/ 



we immediately obtain the following result. 

Lemma 5.6. For any J > 1 and even Nj > 1 (forj = 1,2 ... , J): 



I 



1 



Et E -E 



j = l ^ n<Nj or n^Sj g 
n even 



-i2je]V| V3/ 
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Define for e, c > 

Tj''^ = m) : hj{p"^) > — ^ p prime, m G Z witli m > l| 

Now Tj''^ is finite wlienever < e < 1 (let p*" — > oo in the definition). 
Put 

Mj^e = niax hj{n)rf\ 



note that Mj^e > 1- 



Lemma 5.7. Lei 0<e<l. If n E Sj^e then for m > 1 
in particular, Sj^e is finite. 

Proof. Let M = Mj^^. Suppose that n E Sj^e, and write n = ni ■ n2 
with gcd(ni,n2) = 1, such that if p"^ \\ ni then {p,m) is in Tj'^^ and if 
p"^ II n2 then {p,m) ^ Tj'*^. 
If n2 > 1 then 

hj{n)rf = hj{ni)n\ ■ hj{n2)n2 < M ■ = 1, 

contradicting our assumption n E Sj^e- Thus n2 = 1 and the lemma is 
proved. □ 

This means that the inner sums in Lemma (15. 6p are finite. 

Algorithm 5.8. This results in the following method for computing a 
lower bound for /5. 

(1) Choose J > 1, and perform the following three steps for j = 
1,2,..., J. 

(2) Choose an even integer Nj > 1 large enough. 

(3) Compute 4 J2"<^] l^'jiji); for example from the definitions ( (5. 5\) 

even 

and [5^) . 

(4) Determine the set Sj^e follows: 

(4a) Choose e = cj and determine Tj^e,i- Then compute Mj^e = 
max„ hj{n)n'^ by choosing the product of the worst value for 
{p, m) E Tj^e,i for each prime p occurring in this set. 

(4b) Choose c = Cj and determine Tj^e,c- 

(4c) Determine Sj^e, that is, the positive integers n built up from 
prime powers p"^ with {p,rn) in Tj^e,c, for which hj{n) > 
l/r^^ 

and compute i "ss^,^ 
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(5) Compute the error term (2jeiVJ) ^ ■ y^j . 

(6) Take the sum of results from Steps 3 and 4, o-nd subtract the 
sum of the results of Step 5, taken over j = 1, 2, . . . , J. 



j 


e 


#S 


main 


contribution 


S 


error 


1 


1 





0.508058 






4.18 - 


10- 


12 


2 


0.75 


71678431 


0.134230 


4.7096 ■ 10- 


12 


2.99 


- 10- 


-9 


3 


0.60 


139189128 


0.048944 


8.949 ■ 10" 


12 


3.276 


- 10" 


-7 


4 


0.48 


93183633 


0.020684 


9.7488 ■ 10- 


12 


2.462 


- 10- 


-6 


5 


0.35 


10201152 


0.009564 


-9.1679- 10- 


12 


2.66 


- 10- 


-5 


6 


0.28 


27662520 


0.004706 


-1.95-10- 


12 


7.89 


- 10- 


-5 


7 


0.20 


24415897 


0.002425 


-3.315-10- 


12 


3.31 


- 10- 


~4 


8 


0.15 


65291514 


0.001295 


5.907-10- 


12 


7.26 


- 10- 


-4 


9 


0.03 


7466778 


0.000711 


-9.511 - 10- 


12 


2.59 


- 10- 


-2 


2^j=i 






0.729906 


< 10- 


lu 


1.1625 


- 10- 


-4 



Example 5.9. Take = 10^ The table lists for j = 1, 2, . . . , 9 the 

values of e, j^S, M, the main term from Step (3) in the algorithm, the 
contribution from S in Step (4) and the error term from Step (5). 

As an indication of the size of the numbers involved: the largest 
element of 5*2,0.75 is the 24-digit product of the 18 primes less than 62. 
Since the error term for j = 9 exceeds the contribution of the main 
term, we have not included this one in the final sum. 

As a consequence, we find that P > 0.728743. 

Corollary 5.10. 

A < -0.030. 
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